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1.  OBJECTIVES 

The  principal  objective  of  this  project  is  to  create  an  X-Windows-based  graphics  tool 
to  compute  rapidly  and  efficiently,  synthetic  seismograms  for  laterally  heterogeneous, 
two-dimensional,  isotropic  velocity  models  using  the  Gaussian  beam  method.  The  existing 
Gaussian  beam  software  is  written  in  Fortran  code  and  can  be  very  labor  intensive  to  use. 
Our  goal  is  to  construct  an  X-Windows  Graphical  User  Interface  (GUI)  which  will  elim¬ 
inate  much  of  the  tedium  of  introducing  lateral  heterogeneity  into  two-dimensional  velo¬ 
city  models. 

In  this  report,  we  describe  the  overall  architecture  of  the  modules  and  how  they  inter¬ 
face  with  software  already  developed  under  the  rabric  of  the  NMRD.  A  brief  explanation 
of  how  Gaussian  beam  seismograms  are  computed  is  included  for  clarity.  Finally  we  con¬ 
clude  with  an  outline  of  progress  to  date  and  our  objectives  for  the  remainder  of  the  con¬ 
tract. 

2.  FUNCTIONAL  OUTLINE 

The  functional  flow  for  computing  Gaussian  beam  seismograms  and/or  calculating 
traveltimes  through  heterogeneous  media  is  shown  in  Figure  1.  The  first  step  is  either  to 
create  an  input  model  from  scratch  or  to  access  a  fully  two-dimensional  model  which  has 
been  created  previously.  The  former  is  generally  done  by  beginning  with  a  one¬ 
dimensional  model  and  extending  it  into  a  second  dimension.  Once  in  this  extended  form 
the  user  may  impose  an  overall  heterogeneity,  such  as  ellipticity  in  the  case  of  a  global 
model,  or  a  localized  heterogeneity,  such  as  a  sedimentary  basin  with  low- velocity  lens  in 
the  case  of  a  regional  model.  Additional  heterogeneity  may  be  created  by  manually  mani¬ 
pulating  the  model  elements  with  a  mouse.  Whether  created  or  read  in  and  modified,  the 
model  may  be  stored  at  the  end  of  this  step. 

The  velocity  model  is  specified  as  a  series  of  knotpoints  and  triangles.  Each  knot- 
point  fixes  Vp,  Vj,  and  p  at  a  point  in  space.  Because  the  velocity  gradient  is  assumed  to 
be  linear  between  each  knotpoint,  the  velocity  is  effectively  specified  fully  in  two  dimen¬ 
sions.  (One  exception  is  at  discontinuities:  there  two  knotpoints  are  spatially  co-located 
and  specify  the  velocity  and  density  on  each  side  of  the  discontinuity.)  Knotpoints  are 
grouped  into  triplets  to  form  triangles.  A  value  for  and  Q^,  the  P-  and  5-wave 
attenuation,  is  assigned  to  the  space  enclosed  by  each  triangle.  The  program  tracks  which 
triangles  share  knotpoints  and  are  therefore  "neighbors." 

Under  these  linear  gradient  conditions,  there  is  an  analytical  solution  for  the  raypath 
across  a  triangle.  Raytracing  through  the  model  is  accomplished  by  stepwise  tracing 
analytically  through  each  component  triangle  along  the  raypath.  Anelasticity  is  accounted 
for  by  computing  a  t*  operator  using  and  from  the  triangles.  This  is  the  second 
step  shown  in  Figure  1 .  Results  from  this  step  are  also  stored  for  later  reference. 

The  essentials  for  traveltime  calculation  or  seismogram  computation  are  now  com¬ 
plete.  One  should  remember  that  in  the  Gaussian  beam  method,  it  is  not  necessary  to 
compute  rays  which  travel  directly  from  source  to  receiver.  Rather,  it  is  sufficient  to  com¬ 
pute  a  number  of  rays  which  originate  at  the  source  and  terminate  within  several 
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Functional  Analysis 


Figure  1. 


wavelengths  of  the  receiver.  Therefore,  to  obtain  a  traveltime  or  seismogram,  one  must 
first  specify  the  source-receiver  geometry  and  phase,  or  component  phases  in  the  latter 
case. 
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System  Architecture 


Figure  2. 


This  is  represented  by  the  display/query  box  at  the  bottom  of  Figure  1.  Since  the 
prerequisites  have  all  been  computed,  one  may  repeatedly  access  the  model  and  ray  infor¬ 
mation  to  obtain  traveltimes  and  seismograms  for  different  receivers  and  focal  mechan¬ 
isms. 

3.  SYSTEM  ARCHITECTURE 

How  the  functional  capabilities  are  realized  is  outlined  in  Figure  2.  In  this  diagram, 
rectangles  represent  programs  with  an  X-windows  graphical  component,  ellipses  enclose 
the  names  of  background  processes,  and  ISIS  is  the  name  of  the  interprocess  communica¬ 
tion  (IPC)  module  employed  by  the  NMRD  system.  The  modules  being  developed  for 
this  project  are  Xgb  and  GBseis. 

It  is  Xgh  which  is  at  the  heart  of  this  system.  The  X-windows  interface  allows  the 
user  to  create  or  modify  two-dimcn.sional  earth  models.  If  a  velocity  structure  is  to  be 
assembled  from  scratch,  the  user  is  presented  with  a  suite  of  one-dimensional  models. 
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either  global,  regional,  or  customized  by  the  user.  The  one-dimensional  model  is  recast 
into  triangles  extending  laterally.  Alternatively,  the  user  may  read  in  a  model  created  with 
Xgb  on  a  previous  run.  This  is  accessed  from  where  it  is  stored  on  disc  via  the  Oracle 
Relational  Database  Management  System  (RDBMS)  running  at  the  Center  for  Seismic  Stu¬ 
dies.  A  prototype  application-specific  database  table  to  handle  Gaussian  beam  model 
storage  is  shown  in  Table  1  as  relation  gbmodel. 

Once  the  models  are  input,  they  may  be  modified  in  several  ways.  One  way  is  to 
apply  a  spatial  filter  to  all  knotpoints.  Examples  of  situations  for  which  this  may  desirable 
are  to  account  for  ellipticity  in  global  models,  or  for  synclines  or  inclined  layers  in  local 
or  regional  models.  The  outcome  of  the  filtering  process  is  displayed  immediately,  so  the 
user  can  inspect  to  see  if  the  results  are  satisfactory,  and  if  not,  can  reverse  the  process. 

Additionally,  the  user  may  introduce  heterogeneity  via  the  mouse  by  manipulating  the 
spatial  location  of  knotpoints  or  by  selecting  a  triangle  or  series  of  triangles  and  allowing 
the  user  to  modify  the  properties  of  those  selected.  This  includes  the  velocities  and  densi¬ 
ties  assigned  to  the  constituent  knotpoints  as  well  as  the  attenuation.  The  functionality  is 
provided  as  an  overlay  of  the  filtering  mechanism.  That  is  to  say,  mouse  manipulation 
may  be  done  before  or  after  a  filtering  operation. 

Once  the  model  is  adjusted  to  the  satisfaction  of  the  operator,  raytracing  may  be  per¬ 
formed  with  Xgb.  This  is  accomplished  by  specifying  the  initial  ray  parameter  of  the  ray 
and  by  describing  the  discontinuity  interactions,  namely  what  discontinuities  the  ray 
encounters  and  whether  it  is  reflected,  transmitted,  or  converted  there.  The  user  has 
mouse  control  of  the  source  location,  and  every  time  the  source  position  is  changed,  rays 
are  automatically  retraced  through  the  medium.  The  operator  has  the  option  of  modifying 
the  model  and  retracing  the  rays.  Therefore  iteration  of  model  alteration,  source  position 
change,  and  ray  specification  can  be  done  until  the  operator  wishes  to  preserve  the  results. 
Like  the  model,  raytracing  results  are  stored  in  a  file  for  later  reference,  and  a  pointer  to 
this  file  location  is  passed  to  the  Oracle  RDBMS.  A  prototype  table  for  storing  this 
information  is  shown  in  Table  1  as  relation  gbrays. 

The  computational  functions  of  Xgb  are  complete  at  this  stage.  Traveltime  calcula¬ 
tion  and  synthetic  seismogram  computation  are  performed  by  the  non-graphic  server  pro¬ 
cess  GBseis.  This  is  accessed  via  IPC  messages  instructing  it  what  function,  traveltime  or 
synthetic  seismogram,  to  fulfill.  If  the  function  is  traveltime  calculation,  GBseis  will 
return  its  results  via  IPC  reply.  If  it  is  seismogram  computation,  the  results  will  be  written 
to  disk,  registered  in  the  database,  and  an  IPC  acknowledgment  sent.  Xgb  will  be  able  to 
form  the  proper  IPC  messages  to  prompt  GBseis  to  perform  the  described  calculations,  and 
specifications  of  the  IPC  messages  will  be  provided  so  that  other  processes  will  be  able  to 
do  so  as  well. 

By  design,  there  is  no  capability  in  either  Xgb  or  GBseis  to  display  the  resulting  syn¬ 
thetic  seismograms.  In  keeping  with  the  NMRD  goal  of  modular  design  and  distributed 
processing,  this  task  is  left  to  modules  such  as  the  Analyst  Review  Station  (ARS)  or  Geo¬ 
tool,  both  under  NMRD  development. 
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delation: 

Description: 

gbniudcl 

Gaussian  beam  models 

attribute 

name 

field 

no. 

storage 

type 

external 

format 

character 

positions 

attribute 

description 

modelid 

1 

i4 

i8 

1-8 

model  id 

dir 

2 

c64 

a64 

10-73 

directory 

dfile 

3 

c32 

a32 

75-106 

data  file 

commid 

4 

i4 

i8 

108-115 

comment  id 

Iddate 

5 

date 

al7 

117-133 

load  date 

Relation: 

Description: 

gbrays 

Gaussian  beam  raytracing  results 

attribute 

name 

field 

no. 

storage 

type 

external 

format 

character 

positions 

attribute 

description 

modelid 

1 

i4 

i8 

1-8 

model  id 

dir 

2 

c64 

a64 

10-73 

directory 

dfile 

3 

c32 

a32 

75-106 

data  file 

xsrc 

4 

f8 

fl6.3 

108-123 

x  SIC  coordinate 

zsrc 

5 

f8 

fl6.3 

125-140 

z  sre  coordinate 

commid 

6 

i4 

i8 

142-149 

comment  id 

Iddate 

7 

date 

al7 

151-167 

load  date 

Table  I. 


4.  PROGRESS 

The  basic  core  of  the  module  Xgb  is  nearly  complete.  The  code  for  model  construc¬ 
tion  and  kinematic  raytracing  has  been  translated  from  Fortran  to  C  so  that  it  may  be  more 
smoothly  integrated  into  the  X-windows  display  routines.  The  user  may  now  select  from 
a  suite  of  one-dimensional  models,  the  starting  model  that  is  projected  into  a  second 
dimension  and  displayed  using  NMRD  plotting  widgets  developed  elsewhere. 

Such  a  display  is  shown  in  Figure  3.  Here  a  simple  four-layer  velocity  model  has 
been  extended  laterally.  At  this  stage,  all  knotpoints  at  a  given  depth  share  the  same  velo¬ 
city  and  density  as  befits  the  one-dimensional  model  from  which  they  were  generated. 
What  cannot  be  shown  in  this  black-and-white  representation  is  that  there  are  discontinui¬ 
ties  between  the  first  and  second  layers  and  between  the  second  and  third  layers.  On  a 
color  display,  these  discontinuities  are  shown  in  a  different  color  to  indicate  they  are 
different  from  normal  triangle  sides. 
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Figure  4 


r 


6 


Teledyne  Geotech 


F29601-91-C-DB04 


23  August  1991  -  4  April  1992 


Figure  5 


The  capacity  to  warp  the  model  into  a  simple  synclinal  or  anticlinal  form  is  now 
available.  The  results  of  applying  this  simple  filter  may  be  seen  in  Figure  4.  Here  the 
one-dimensional  model  of  Figure  3  has  been  filtered  to  deform  it  into  a  broad  synclinal 
structure.  The  properties  of  the  knotpoints  have  not  been  changed  -  only  their  vertical 
coordinates.  As  stated  above,  this  is  just  one  of  several  filters  the  operator  will  be  able  to 
apply  in  the  completed  version. 

The  code  for  accomplishing  kinematic  raytracing  is  complete.  Figure  5  shows  rays 
reflecting  off  a  discontinuity  of  the  model  shown  in  Figure  4.  The  triangles  are  omitted 
from  the  display  here  for  clarity  —  only  the  rays  and  the  discontinuities  are  shown.  Before 
synthetic  seismograms  can  be  computed,  the  dynamic  aspects  of  the  raytracing  must  be 
incorporated.  This  is  a  straightforward  extension  of  the  raytracing  functions  already 
coded. 

5.  FUTURE  PLANS 

There  are  three  principle  developments  which  must  be  accomplished  in  the  near  temi 
for  this  system  to  function:  (1)  the  dynamic  raytracing  must  be  coded  into  Xgb,  (2)  the 
code  for  the  GBseis  module  must  be  translated  from  Fortran  to  C  and  IPC  hooks  installed, 
and  (3)  the  Xgb  interface  must  be  modified  to  allow  the  user  to  manipulate  individual 
knots  and  triangles.  At  this  writing,  the  logic  for  (3)  had  been  conceived,  and  we  antici¬ 
pate  only  two  or  three  man-days  for  it  to  be  completed.  Likewise,  (1)  involves  another 
Fortran  translation,  but  with  the  C  structures  already  implemented,  this  should  not  require 
more  than  1-2  man-weeks.  The  greatest  task  may  be  (2),  but  we  expect  tliat  the  structure 
of  C  is  so  much  better  suited  for  computations  of  this  kind  that  a  substantial  amount  of 
resources  will  not  be  required  to  complete  that  task. 
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When  these  problems  are  solved,  refinements  to  the  Xgb  will  be  added  to  make  the 
module  easier  to  use.  Included  in  this  are  plans  to  (1)  allow  phases  whose  rays  are  to  be 
computed  to  be  selected  from  a  menu  list,  (2)  specify  source  type  as  either  explosion  or 
double  couple,  and  if  the  latter,  allow  the  user,  via  graphics,  to  alter  the  orientation  of  the 
fault  planes  and  have  the  changes  be  reflected  in  the  computed  seismograms,  and  (3)  allow 
the  user  to  project  "target"  receivers  onto  the  model  and  check  the  ray  density  of  rays  ter¬ 
minating  within  a  fixed  number  of  wavelengths  of  the  target. 
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